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Abstract We present a complementary study to a new model for generating magnetic fields of cosmological in- 
terest. The driving mechanism is the photoionisation process by photons provided by the first luminous sources. 
Investigating the transient regime at the onset of inhomogeneous reionisation, we show that magnetic field am- 
plitudes as high as 2 x 10 -16 Gauss can be obtained within a source lifetime. Photons with energies above the 
ionisation threshold accelerate electrons, inducing magnetic fields outside the Stromgren spheres which surround 
the ionising sources. Thanks to their mean free path, photons with higher energies propagate further and lead to 
magnetic field generation deeper in the neutral medium. We find that soft X-ray photons could contribute to a 
significant premagnetisation of the intergalactic medium at a redshift of z — 15. 
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1. Introduction 

Various observational techniques have been employed to 
detect the presence of magnetic fields in very different 
structures and on very different scales. On the largest 
scales, the measured amplitudes of the magnetic field are 
generally in the range of a few tent hs microgauss to tens 
of microgauss, bo th in galaxies (e . g.lBeck et al.ll l996tl an d 
on cluster scales (Kr onberdl2001t ICarilli fc Tavlon l2002). 
Faraday rotation measurements have permitted to ob- 
tain an upper limit for the intergalactic magnetic field of 
I ^igm | ^5 1 ~ 9 G, assumin g a maximal field reversal scale 
of 1 Mpc teronberdllOol . 

The origin of those magnetic fields remains largely un- 
known. It has been soon realised that cosmic magnetic 
fields could not have been born directly with such high 
amplitudes. Had it been so, the structure formation his- 
tory would have been drastically modified, since a field of 
10~ 9 G on megaparsec scales is strong enough to induce 
nonline ar matter density fluctuations l)Wassermanlll978l 
see also lKim et alJll99fift . Clearly then, cosmic magnetic 
fields must first arise in the form of weak seeds, the am- 
plitude of which is then amplified by some powerful mech- 
anism. 

A plethora of models for the creation of cosmic mag- 
netic fields is available in the literature, but none of the 
mechanisms envisaged seems completely satisfactory. On 
the one hand, there is a large variety of models relying on 
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the physics of the early universe, but the predicted seed 
amplitude is strongly model dependent, and spans conse- 
quently a rather large interval (10~ 65 G to 10 -9 Gauss 
typically - for rece nt 

2nmt lWiHrowll2nn2l l 
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On the other hand, after recombinatio n, most of the 
astro-physical mechanisms pr oposed (e.g. IPudritz fc Silkl 
ll989t lD avies fc Widrowll200(t are based on the Biermann 
battery mechanism l|Biermannlll95(tl . first investigated to 
account for stellar magnetism. Applications of the battery 
mechanism to various cosmological contexts lead to seed 
fields of the order of 10 -20 — 1CU 17 G on scales of a few 
kiloparsecs to a few hundreds of kilo parsecs (for the con- 
text of reionisation in pa rticular, see ISubramanian et alJ 
ll994UOner]in et alJ^OOrf Since those models rely on the 
same mechanism, namely differential acceleration of elec- 
trons (vs. ions) due to thermal pressure, acting on similar 
scales, it is not surprising that the range obtained for the 
initial magnetic field is rather narrow (as compared to 
early universe models). Those fields are usually thought 
to be sufficient to account for present day galactic fields 
provi ded a subsequ e nt a mplification by galactic dynamo. 

In lLanger et~afl ^200.^ . we proposed an efficient mag- 
netogenesis model in which the driving force is the radi- 
ation emanating from the first luminous objects. Under 
the assumption of stationary regime, we found that 
anisotropic and inhomogeneous radiation flux leads, at 
the end of cosmological reionisation, to magnetic fields 
as high as 1CU 12 — 1CU 11 Gauss on protogalactic scales. 
Moreover, we also found that magnetic fields on small 



scales are strongly clamped. Ihe advantage ol our model 
resides mainly in two facts. First, the driving force is the 
radiation drag which provides an acceleration of charged 
particles inversely proportional to their mass cubed. As 
a consequence, electric currents develop millions of times 
more efficiently as compared with Biermann battery mod- 
els which rely on thermal pressure. Second, the properties 
of the magnetic field power spectrum are directly related 
to the shape of the power spectrum of the matter den- 
sity fluctuations. This guarantees that large scale fields 
are favoured relative to small scale fields. 

In the present article, we extend our first model to the 
case of the transient regime. More specifically, we focus 
on the generation of magnetic fields at the onset of reion- 
isation. The driving process for charge separation and ac- 
celeration is photoionisation itself. The formalism and the 
results are given in Sects. |21 and [3J Discussion and conclu- 
sions are presented in Sects. 0] and [SJ 

2. Formalism 

We study the generation of magnetic fields due to electric 
currents induced by the acceleration provided by radia- 
tion drag. The relevant equations are the usual Maxwell 
equations 



(see lFeebles et al.ll/UUUl atter lbeaeer et al.1120 00). and the 
plasma frequency is uj p ~ 26rads _1 , which is much larger 
than the electron ion collision rate, v c ~ 1.35 x 10 _5 s _1 
(at a mean gas temperature T ~ 0.85 K, as redshiftcd 
from T ~ 4000 K at z = 1100 to z = 15). Finally, esti- 
mating the photon flux at the hydrogen ionisation thresh- 
old as cj> ~ 4.56 x 10 31 cm~ 2 s~ 1 (at a distance of 10 2 
kpc of a 10 12 L© source, see eq. ijlUl below), we obtain 



3 x 10 7 



considering photon-electron Thomson 



scattering. 

Thus, the development of electric currents from j = 
can be decomposed in three distinct time sequences. 



First, for times shorter than 



only the source term is 



relevant. Considering it is constant in time, we obtain 
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So initially, the magnetic field grows quite quickly, but this 
regime is valid on a very short time only, i.e. for t < <j~ . 

However, very soon t reaches w" 1 , and in that case 
we need to consider the electric part of Lorentz force, in 
which case the system to solve is 



V • E = d t E = cV x B - Anj 

V ■ B = d t B = -cV x E 



(la) 
(lb) 



where the electric current is j = n e q e (v e — v p ), n e be- 
ing the electron density, q e the electron charge, u e (p) * ne 
electron (proton) velocity, and we have supposed that the 
medium is globally neutral. We also consider the gener- 
alised Ohm's law, 



d t j + (u • V) j = 
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u x B 
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(2) 



where uj p is the electron plasma frequency, v c is the 
electron-ion collision frequency, m e the electron mass, and 
it is the peculiar velocity field of the center of mass of the 
plasma. We included the radiation drag effects as a source 
term, 
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hv 



m e c 



(3) 



defining the photon plasma interaction frequency by v~ t = 
a-y<p , where hv is the photon energy, is the photon 
flux and c 7 the photon-electron interaction cross-section. 
The source current / obviously depends on the ionisa- 
tion state of the medium (see Sect. Initially, all but 
the source term are zero, and we may explore the gen- 
eration of the magnetic fields by considering linear terms 
only. We defer to Sect.E]the discussion of nonlinear terms. 
Equation J3J shows three characteristic frequencies, ui p ~ 



5.64 x 10 4 ne / "rads~ 1 , v c - 50n e r- 3 / 2 s- 1 and u y . For il- 
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lustration, at redshift 15, when the universe is almost neu- 
tral, we have a residual ionisation of n e ~ 2.1 x 10~ 7 cm -3 



d t j = v^I+^-E 

d t E = -Anj + cV X B 
d t B = -cV x E. 

Using Fourier and Laplace transforms, we obtain 



E k = -4nu 



1 — cos [ui T t] 
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3k = 



B k 



t^ p sin (u) T t) 



sin (oj T t) 



k x J fe 



(5a) 

(5b) 
(5c) 

(6a) 
(6b) 
(6c) 



where oj t — (w 2 + c 2 fc 2 ) 
uj p /c ~ 1.82 x 10 



1/2 



ujp for scales k <?C k p = 

Taking for instance the 
residual ionisation after recombination at a redshift of z ~ 
15, we have n e ~ 2.1 x 10 _7 cm~ 3 , and the corresponding 
plasma scale is fc" 1 ~ 1.2 x 10 7 m. After reionisation is 
complete, say at z ~ 6, we have n e ~ 8.1 x 10 -5 cm -3 , and 
the plasma scale is even smaller. Therefore, for all scales 
of interest, the characteristic pulsation is oj t = ui p . 

Because of the intrinsic velocity distribution of the 
electrons, the oscillations shown in eqs. © will vanish by 
Landau damping, essentially within a few t p ~ w" 1 . This 
means that after a few plasma time scales, the electric field 
can be considered constant, the electric current is vanish- 
ing (except perhaps on very small scales, i.e. for k S> k p , 
where j continues to grow linearly) , and the magnetic field 
is increasing linearly with time. 

Finally, the next time scale involved is the magnetic 
field diffusion time, Tjy = Loik -2 / (c 2 v c ), until which the 



evolution of is will be given (from eqs.|Dj) by 

B = 4^tV x/ (7) 

on scales bigger than fc" 1 . 
3. Magnetic field seeds 

We can now estimate the order of magnitude reached by 
the magnetic field in a source lifetime t§. However, in the 
vicinity of an ionising source, the actual physical process 
leading to electron acceleration will depend on the ionisa- 
tion state of the medium. 



3.1. The Stromgren sphere 

The extent of the ionised region around a luminous source 
depends on the recombination rate. If recombinations are 
negligible, the radius of the ionised sphere depends on the 
age tg of the source, 



rs 



3N ph t s 
47r(n H ) 



1/3 



(8) 



where N p h is the source ionising photons emission rate and 
(tih) is the mean hydrogen density. In the opposite case, it 
is determined by the balance between recombination and 
emission of ionising photons, and its value is given by 



37V, 



1/3 



rs = 



ph 



47ra rec C(n H } 2 



(9) 



where a roc = 2.6 x 10 13 cm 3 s 1 is the radiative re- 
combination coefficient a t a gas temperature of 10 4 K 
^Hummer fc Seatoul M flfifll and C = (n|)/(n H > 2 > 1 is 
the clumping factor. 

The recombina tion time scale can be defined as (e.g. 
iMarfaii et alF'MHI. 



[1.17(n H )a r o C C] 



1.26 x lO 8 ^ 1 



0.024 



16 



-3 



yrs (10) 



Depending on its definition and on the model of struc- 
ture fo rmation, very differe nt values for C are ob- 
tained. iBenson et alJ l)200l[) for instance obtain val- 
ues sligh tly above 3 at z = 15, whereas calcula- 
tions by Valagcas & Silk (1999) yield C ~ 40 at the 
same redshift. Also, results from simulations are quite 
sensitive to numerical resolutio n a p oint early em- 
phasized bv iGn edin k. Ostrikerl ljl997|) . For instance, 
ISpringel fc Hernauistl l|2003j) obtain C in the range 2.5 — 6 
at z = 15, with higher values corresponding to higher res- 
olutions. Facing those uncertainties, we adopt C = 10 as 
a fiducial value at z = 15. 

Let us consider quasars as the sources of the ionis- 
ing radiation. The typical lifetime of quasars is not a 



tightly constrained quantity, and values ranging from 10 
to 10 8 years are usually quoted (e.g. lMartinill2004]) . It is 
nevertheless reasonable to assume that quasars spend an 
Eddington time scale, 



t E = 4 x 10 s e yrs, 



(11) 



in their shining phase, and thus live at least that long 
(e ~ 0.1 is the radiative accretion efficiency). Then, com- 
paring eqs. (TTU1) and ifTTjl. we see that for z = 15, the 
recombination time scale is shorter than the quasar life- 
time as soon as C > 3.15. Thus, the radius of the ionised 
bubble can be considered constant. The total time spent 
by the bubble in its expanding phase is simply given by 
the ratio of the recombination time scale to the source life- 
time. Considering ts as the minimum lifetime, this ratio 
is at most of order 



tree 



= 0.315 e^C 



0.024 



l + z 
16 



(12) 



and is less than 13% of the source lifetime at z = 15 for 
t$ ~ 10 8 yrs and C = 10. In the following, we therefore 
adopt the definition of rs given by eq. 

Whichever the photon source, the most relevant pho- 
tons are those that are roughly at the hydrogen ionisation 
threshold. Indeed, for a stellar source, the UV spectrum 
is approximately flat with a cutoff at only a few times the 
Lyman limit frequency. Similarly, quasar spectra invari- 
ably exhibit the so-called Blue Bump at A ~ HOOA (e.g. 
iTelfer et al.ll2002t IShang et al.ll2005[) . Let us then count 
as ionising photons those with energies in the range 13.6 
eV - 29.3 eV (ionisation by photons with higher energies 
is more than ten times less likely to happen) . The ionising 
photons emission rate can be then estimated as 



N, 



ph 



L Eo L f£ 3 <l>{r,E)dE 



j e 
E 



En f^<t>(r,E)dE 



(13) 



where Eq — 13.6 eV, L is the source luminosity and E c is 
the energy at the spectrum cut-off. For quasars, we have 



(r,E) = Mr)(E/E ) a exp(E/E c ), 



(14) 



with a ~ -1.7 an d E c ~ 2 keV l|Telfer et all l2002t 
ISazonov et a"fll2004|) . and 



ph 



0.4- 



L 

En 



7 x 10 55 Li2 s" 1 



(15) 



with L12 = L/10 12 Lq. Notice that all z > 6.1 quasars 
observed to date are extremely luminous and their 7V p h 
estimated from observations l ie in the range [0.6 — 2.2] x 
10 57 s" 1 (see, e.g., Table 1 of lYu fc LiJEooi . Finally, we 
obtain a Stromgren radius of 



r s ~ 683 C. 



-1/3 r 1/3 



10 



v 12 



il b h \ fl + Z 



0.024 



16 



kpc, (16) 



where C w = C/10. 



6.Z. Ionised case 

Within the Stromgren sphere, the plasma is basically fully 
ionised, and the dominant mechanism to accelerate elec- 
trons is simply the Thomson scattering of photons off free 
electrons. Therefore, u 1 = <tt = 6.65 x 10~ 
source current I in this case is defined by 



/ = 



E 



m e c 



- q e n e c. 



Then, directly from eq. J7J), we get 



B 



R mt q e 



ts 



(17) 



(18) 



where R Iot is the scale of inhomogeneities in the ionising 
flux <j). At the ionisation energy threshold, we can quickly 
estimate the relevant photon flux at a distance r from the 
source by taking 



E4>(r,E) = 



L 



47rr 2 



at E = 13.6 eV. The amplitude of the field is then 



B ~ 4.56 x 10 



-28 



-12 



R\ ot rl 10 s yrs 



Gauss 



(19) 



(20) 



where R 2 ot 



i? rot /10 2 kpc and r 3 = r/10 3 kpc. The 
field strength we obtain here is extremely small, and has 
presumably little relevance for seeding the protogalactic 
plasma (see also Sect. 13.411 . Notice that for small r the 
amplitude of the field in eq. i12()ll could be much higher. 
However, at short distances from the source, the plasma 
is likely to be turbulent, and outflows perhaps driven by 
stellar winds, or jets powered by central engine activity, 
may occur. In those situations, pure plasma instabilities, 
such as the Weibel instability, are likely to be much more 
efficient in producing magnetic fields than the model pre- 
sented here. 

3.3. Neutral case 

Beyond the Stromgren radius, the situation is quite differ- 
ent. Indeed, the medium is opaque to UV photons at the 
ionising energy threshold, and the medium is therefore ba- 
sically neutral. However, photons with higher energies do 
propagate further deep into the neutral region and even- 
tually ionise atoms at a distance which depends on the 
ionisation mean free path. In this case, the mechanism of 
electron acceleration is therefore the ionisation process it- 
self. Thus, fi 7 = Cion is now the ionisation cross-section, 
which depends on the photon energy as 



00 



-3 



(21) 



with (j o ~ 7.8 x 10 18 cm 2 , and the source current is now 
defined by 



fn 



E 

' m e 6 



2<lenilC, 



(22) 



E (eV) 

Figure 1. F(E) = f mt (E / E a )~ 3 - 7 exp(E/ E c ) as a func- 
tion of the incident photon energy. 

where «h is the density of neutral hydrogen, and / m t is the 
fraction of the incident photon momentum transferred to 
the photoelectron (notice that we neglected terms of the 
order of m e /m p ). The amplitude of the seed field is now 
estimated to be 



■is- 



(23) 



n e R mt q e 

Using eq. (|14J) as a definition for the ionising photon flux, 
we normalise the spectrum at a distance r by estimating 
the power emitted by the source, i.e. 

' <t>{r)dE = JL, (24) 



47rr 2 



which gives 



<j>o(r) ~ 5.4 x 10" 



L 



E r 2 



(25) 



and in this case, the field amplitude obtained is 
(n H /« e ) L u t s 



B - 7.6 xlO~ 16 .F(£;)- 



10 4 i?|°V| 10 s yrs 



Gauss, (26) 



where F(E) = f mt (E/E a ) 3 ' 7 exp (E/E c ). Notice that in 
the equation above, as compared to eq. I|20|l . r 3 has been 
replaced by rg/lMpc, the Stromgren radius normalised 
to 1 Mpc, since we are interested in the effects of photons 
penetrating into the quasi-neutral medium outside the Hn 
bubbles. Using the definition of the Stromgren radius in 
the amplitude of the field obtained is finally 



eq 



en, 



B 



1.63 xlO- 15 ^)^ 3 ^^ 



L 



1/3 



ts 



R 2 ot 10 8 yrs 



n b h 2 

0.024 



4/3 



1 + Z 

16 



Gauss. 



(27) 



Now, the fraction of momentum transferred to a pho- 
toele ctron depends on the energy of the incident photon 
(e.g. lMassacrieilll996l) , 



fmt 



8 E-E t 
5 E~ 



(28) 



E (eV) 


17.3 


25 


30 


53 


136 


Zmfp (10* P C) 


1.03 


3.12 


5.39 


29.7 


502 


B (10"" G) 


23.1 


12.7 


7.77 


1.3 


0.05 



Table 1. Photon mean free path in the neutral medium as 
a function of their energy and corresponding values of the 
generated magnetic field at z = 15 (assuming a fluctuation 
scale of 100 kpc and a clumping factor of C = 10). 



such that / mt ~ at the ionisation threshold, and the am- 
plitude of the field obviously vanishes. However, as shown 
on Fig. where we plot F(E) as a function of E, as soon 
as we are slightly above the threshold, B > 0, and the 
most favorable photons for magnetic field generation have 
E ~ 17.3 eV. The photon mean free path outside the 
Stromgren sphere is Z m f p — (fTi on ?iH) _1 , i-e. 

and is increasing with energy. This implies that the 
strongest fields are created just within the thin skin of 
the ionised bubble, whereas slightly weaker fields are gen- 
erated on thicker depths. Assuming a clumping factor of 
C = 10, Fig. [21 shows the amplitude of the magnetic field 
we obtain (solid line), as well as the photon ionisation 
mean free path (dashed line), as functions of the photon 
energy. Table ^ summarises a few typical values of the 
magnetic field and of the depth on which it is created. 
As can be seen, the amplitude of the field remains above 
10~ 17 Gauss on kiloparsec scales and decreases down to 
5 x 10~ 19 Gauss at 50 kpc. 

These amplitudes depend on the scale R mt of the inho- 
mogeneities in the neutral gas surrounding the Stromgren 
sphere. The maximum scale that we can consider is there- 
fore of order rg/2 ~ 340 kpc. On the other hand, from 
structure formation, we expect fluctuations of matter den- 
sity on all scales. However, neutral hydrogen fluctuations 
will be present only in haloes in which the virial tempera- 
ture is above 10 4 K. The como ving number density of such 
haloes is roughly 20/i 3 Mpc" 3 l|Mo k Whitell2002t) . which 
corresponds to a proper mean separation of ~ 32 kpc at 
z = 15. Below that temperature threshold, in absence of 
heavy elements, the primordial gas is unable to cool and 
condense. Thus, taking R rot — 100 kpc for illustrative pur- 
poses in eq. (|27|l is fully justified, and the actual value of 
B may differ at most by a factor of ~ 3. 

In addition, the amplitude of the field we obtain de- 
pends on the actual ratio nn/n e at redshift z = 15. 
This ratio is essentially given by the residual ionisation 
state of the universe from recombination era. By taking 
n e /riH ~ 10 -4 in eq. I|27|l . we may have slightly overesti- 
mate d the recombination efficiency. In fact, recent calcula- 
tions l|Seager et al.l200oT) including non-e q uilibrium effect s 
give rather n e /riH ~ a fewlO -4 (see also lHu et al]ll995)) . 

Finally, note also that the amplitude of the field is 
only weakly dependent on the luminosity of the source, 
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Photon energy (eV) 

Figure 2. Amplitude of the magnetic field (solid line) ob- 
tained and photon mean free path (dashed line) as func- 
tions of the incident photon energy. 

meaning that quasars with L = 10 10 — 10 13 Lq will pro- 
duce fields of similar strength. This broadens the applica- 
bility of the mechanism presented here, in that it is not 
restricted to the most powerful, and rare, sources. 

3.4. Source term 

Examining eq. , it is obvious that the magnetic field will 
be non-zero only if the source term has a rotational part. 
This is actually the case if the medium in which photons 
propagate is inhomogeneous. Indeed, if we decompose the 
radiation flux into a mean background <fio and an inhomo- 
geneous component Sip, 

4> = 0o + Sip = (1 + /O0o, (30) 

we can relate the fluctuations of the flux to the fluctua- 
tions of the optical depth t; according to 

fx = e~ Tl - 1 (31) 

and thus also to the matter density inhomogeneities 
present in the plasma, since 

ti = a 7 n e / 5(r)dr n (32) 
Jo 

where n e is the mean electron number density, I is the 
photon propagation path, and S(r) — (n e (r) — n e )/n e is 
the density contrast. Finally, in the optically thin limit, 

/; ~ -G 1 h e I S(r)dr n , (33) 
Jo 

and thus 

Vx/« (V/, x 0o + fiV x O ) ■ (34) 

Obviously, since 4>q depends only on \r\, the second term 
above vanishes. The first term however, has no reason to 
be zero, since at a given distance I outside the Stromgren 
sphere, the neutral medium is inhomogeneous, and the 



optical depth has a non trivial angular dependence. Inside 
the ionised bubble, on the contrary, the ionisation process 
has probably washed out most of the density fluctuations. 
These effects add to the tiny amplitude obtained in eq. 
H2Q(I , meaning that radiation drag most probably does not 
create significant fields within the Stromgren sphere. 

4. Discussion 



In lLanger et"afl {2003), we presented a first mechanism 
for the generation of magnetic fields. In that paper, we 
explored, in the steady state, the effects of radiation drag 
provided by ionising sources at the end of reionisation, i.e. 
when the inhomogeneous universe is fully ionised. The am- 
plitudes that could be reached in that case are as high as 
10~ 12 Gauss on protogalactic scales. Here, we considered 
the early stages of inhomogeneous reionisation. The main 
difference in this situation is that the dominant interaction 
between electrons and photons is photoionisation, of which 
the cross-section is orders of magnitude higher, in a suit- 
able photon energy range, than that of Thomson scatter- 
ing (which dominates in the fully ionised case). Therefore, 
we might have expected to get even higher magnetic field 
amplitudes than in the ionised case, albeit the momentum 
fraction transferred from photons to electrons is smaller 
than one for photon energies below 36.3 eV. However, the 
amplitudes we obtain here are in fact not higher than in 
the fully ionised case. This is mainly due to the combined 
effects of a finite source lifetime and of the early dynamical 
importance of nonlinear terms. 

For initially vanishing fields, nonlinear effects in eq. 
(J2J will not alter the generation scenario outlined above. 
They will, however, play a major role for the field ampli- 
fication once a critical value is reached. Considering the 
fluid velocity field due to structure formation, it is easy 
to show that the first nonlinear term to enter into play in 
the generalised Ohm's law is u x B. 

In the linear regime, the structure formation process 
creates a velocity field 



2 f(tt) 

u = r^-j^g with V • g = inGpS 



(35) 



with /(f2) ~ fi 0,6 . This term enters into play in the usual 
induction equation for the magnetic field as 



d t B = V x (u x B) + AttAt V x I 



(36) 



and dominates over the source term as soon as B ~ 
47r(c/u)(^ 7 /u; 2 )/, i.e. at a time t n \ ~ R vot /u (from eq. 
0) . On a given scale, the velocity field is best described by 
its Fourier mode 



\ u k\ = ; Ok, 



(37) 



and for k ~ l/i? rot , 

t nl ~ (Hf^Sk)- 1 ~ 4.1 x 10 s S^ 1 yrs (38) 



at z = 15 with 11 = 1, \l m = U.27, and h = (J. 71. 1ms 
shows that nonlinear terms enter into play basically im- 
mediately after the period of magnetic field generation. 
It implies that magnetic field amplification by shear and 
adiabatic contraction sets-in directly from the seeds cal- 
culated in Sect. 13.31 without any decrease in amplitude by 
diffusion or dilution due to the expansion of the universe. 

Finally, as seen from eq. I|29|) photons with higher en- 
ergies propagate further, so it would be tempting to calcu- 
late the magnetic field produced by photons with E ~ 400 
eV for instance. Notice however that photons with an en- 
ergy E > 350 eV travel freely distances longer than 880 
kpc, which corresponds to the mean proper distance be- 
tween those rare sources which are embedde d in 3 .9a dark 

l l2(102h . Those 



matter haloes at z = 15 (e.g. iMo fc Whit 
photons therefore do not contribute to the generation of 
magnetic fields on such distances, as there are no more 
privileged directions along which electric fields could be 
generated, but they rather contribute to establish a uni- 
form soft X-ray back ground which partially i onises the 
intergalactic mediu m l|Venkatesan et all l200ll lohl l200ll 
iMadau et alll2004|) . Nevertheless, photons with energies 
~ 280 eV, travelling roughly half the mean source sepa- 
ration, could lead to the creation of magnetic fields with 
amplitudes of the order of 4 x 10 -20 Gauss, premagnetising 
therefore the intergalactic medium at a redshift of z = 15. 

5. Conclusion 

We reported here on a new model for the generation of cos- 
mological magnetic fields. We extended the mechanism of 
charge acceleration by radiation drag, which we first ex- 
plored in a previous a rticle in the station ary regime, at 
the end of reionisation ijLanger et alJl2003|) . In the present 
paper, we considered the transient regime in the case of 
a neutral, inhomogeneous universe at the onset of reioni- 
sation, at a redshift of z = 15. The driving process is the 
photoionisation itself which leads to the creation of mag- 
netic fields around the Stromgren spheres produced by the 
ionising sources. In a source lifetime, the magnetic field 
amplitudes we obtain are in the range 2 x 10~ 16 — 4 x 10~ 20 
Gauss on 100 parsecs to 880 kiloparsecs, i.e. roughly from 
the Stromgren sphere skin to the mean distance between 
sources at z = 15. The growth of the amplitude is limited 
by the source lifetime is, but on the other hand, nonlinear 
terms, driven by velocity fields due to structure formation, 
become of significant importance at about ts- This means 
that the fields begin their amplification by shear and adi- 
abatic contraction directly, without decrease in amplitude 
due to diffusion or by the expansion of the universe. 
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